Voter models on weighted networks 

Andrea Baronchelli, 1 Claudio Castellano, 2,3 and Romualdo Pastor-Satorras 1 

1 Departament de Fisica i Enginyeria Nuclear, 
Universitat Politecnica de Catalunya, 
Campus Nord E>4, 08034 Barcelona, Spain 
2 Istituto dei Sistemi Complessi (ISC-CNR), 
Via dei Taurini 19, 1-00185 Roma, Italy 
3 Dipartimento di Fisica, "Sapienza" Universita di Roma, 
P.le A. Moro 2, 1-00185 Roma, Italy 
(Dated: June 30, 2011) 

Abstract 

We study the dynamics of the voter and Moran processes running on top of complex network 
substrates where each edge has a weight depending on the degree of the nodes it connects. For 
each elementary dynamical step the first node is chosen at random and the second is selected with 
probability proportional to the weight of the connecting edge. We present a heterogeneous mean- 
field approach allowing to identify conservation laws and to calculate exit probabilities along with 
consensus times. In the specific case when the weight is given by the product of nodes' degree raised 
to a power 6, we derive a rich phase-diagram, with the consensus time exhibiting various scaling 
laws depending on 6 and on the exponent of the degree distribution 7. Numerical simulations give 
very good agreement for small values of \6\. An additional analytical treatment (heterogeneous 
pair approximation) improves the agreement with numerics, but the theoretical understanding of 
the behavior in the limit of large \6\ remains an open challenge. 

PACS numbers: 89.75.-k, 89.75.Hc, 05.65. +b 
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I. INTRODUCTION 



Many technological, biological and social networks are intrinsically weighted. Each link 
has associated an additional variable, called weight, which gauges the intensity or traffic of 
that connection, and that can exhibit widely varying fluctuations EUjJ. The presence of 
weights is extremely relevant in some scenarios (e.g. in the case of transport in a network in 
which weights measure bandwidth or capacity), and it must therefore be taken into account 
explicitly. Some results have already been produced in this direction, dealing, among other 

5], epidemic spreading [(3, B|, general equilibrium and 
{J , or glassy dynamics |10|. Here we present a detailed 



problems, with diffusive processes 
non-equilibrium phase transitions [8 



investigation of the ordering dynamics of voter-like models on weighted networks 11). 



The voter model 



12 



13| and the Moran process [l^] are simple examples of ordering 



dynamics, which allow to understand how natural systems with an initial disordered config- 
uration are able to achieve order via local pairwise interactions. Both models are described 
in terms of a collection of individuals, each endowed with a binary variable Sj, taking the 
values ±1. The elementary step consists in randomly choosing a first individual and then 
(again randomly) one of her nearest neighbors. In the voter model the first individual will 
copy the state of her neighbor. In the Moran process, on the other hand, she will transmit 
her own state to the neighboring node, which will adopt it. In both cases, starting from a 
disordered initial state, the iteration of the elementary step leads to the growth of correlated 
domains and, in finite systems, to an absorbing uniform state in which all individuals share 
the same state (the so-called consensus). In a social science context, the voter model repre- 
sents thus the simplest model of opinion formation in a population, in which individuals can 
change their opinion as a function of the state of their nearest neighbors In the same 
way, in a biological context, the Moran process represents the elementary example of two 
species competing (through reproduction and neutral selection) for the same environment 
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The voter and Moran processes are equivalent on regular lattices and on the complete 
graph, but if the pattern of connections is given by a complex (unweighted) topology they 
behave differently, since the order in which interacting individuals are selected becomes rel- 



evant 



16 



171 ]. Moreover, the time to reach consensus scales with the system size in different 



ways depending on the degree distribution of the network 



18 W 

Considering a weighted 
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topological substrate adds, as we will see, a richer and more complex phenomenology. Addi- 
tionally, the case of weighted networks allows to model very natural settings. In the context 
of social sciences, for example, weights can reflect the obvious fact that the opinion of a 
given individual can be more easily influenced by a close friend rather than by a casual 
acquaintance. On the other hand, in an evolutionary scenario, weights allow to gauge the 
effects of heterogeneous replacement rates in different species. 

On weighted networks, at each time step a vertex % is selected randomly with uniform 
probability; then one nearest neighbor of i, namely j, is chosen with a probability propor- 
tional to the weight Wij > of the edge joining i and j. That is, the probability of choosing 
the neighbor j is 

p.. = Wi i (1) 

Vertices i and j are then updated according to the rules of the respective models. With this 
definition, the models considered represent the natural extension for ordering dynamics on 
weighted networks (and in particular of the voter model) of the generalized Moran process 



proposed in Refs. 



21 



22j , in which dynamics was defined as a function of a set of arbitrary 



interaction probabilities Py. In our case, however, the fact that these interaction proba- 
bilities arise from the normalized weights arriving at a vertex imposes some restrictions to 
the possible values of Py and yields therefore different outcomes and interpretations. Also, 



it is worth noting three recent publications 
models. 
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25|, dealing with related, but not identical, 



Adopting the heterogeneous mean-field (HMF) approximation 26|, |27|] we will assume 
that the weight between vertices i and j depends only on the degrees at the edges end- 
points, namely fc, and kj, and therefore we can write = g(ki, kj)atj, where a^- is the 
adjacency matrix and g(k, k') is a positive definite, symmetric function. The application of 



HMF theory and the backwards Fokker-Planck formalism 



I8M20I]. allows us to derive ana- 



lytical expressions in degree uncorrelated networks for the main relevant quantities (namely 



exit probability and consensus time 



ll|) in a more transparent way than in Refs. 2l|, [22] 



In order to allow for closed mathematical solutions of the models, we will further specify the 
function g to be given by the product of two independent functions g(k, k') =_Qs{k)g s {k'), an 
assumption motivated by empirical observations in real weighted networks Specializing 
both models to the case of networks with power-law distributed degrees and edge weights 
given by multiplicative powers the endpoint degrees, g s (k) = k e , a very rich phase-diagram 



is obtained, with several different scaling regions of the consensus time with the network size 
N. A numerical check of the analytical predictions reveals a good agreement in some re- 
gions of the parameters space and noticeable discrepancies in others. In order to gain insights 
into the observed numerical disagreement, we apply an improved mean-field approach, the 



heterogeneous pair approximation |28|, which turns out to provide better agreement with 
numerics for small 9 but is still not able to solve the problems for large 9. The qualita- 
tively different nature of the dynamics for large 9 is briefly discussed and its understanding 
identified as an intriguing challenge for future work. 

II. HETEROGENEOUS MEAN-FIELD THEORY 



In this Section, we perform a theoretical analysis of the voter and Moran processes on 
weighted networks within a HMF approximatio n |26 |. extending the Fokker-Planck formal- 
ism developed for the unweighted case in Refs. jl8l-l20|. Let us consider the models defined 
by the interaction probability Eq. (TTJ), where the network weights take the form 



(2) 



ited networks is to follow 
5] . The key point consists 



The simplest way to extend the Fokker-Planck approach to weig 
the annealed weighted network approximation introduced in Ref. 
in considering the degree coarse-grained interaction probability P w (k — > k'), defined as the 
probability that a vertex of degree k interacts with a nearest neighbor vertex of degree k! . In 
unweighted networks, this probability simply takes the form of the conditional probability 
P(k'\k) that a vertex of degree k is connected to a vertex of degree k' 29]. In networks with 
weights given by Eq. ([2]), the interaction probability of the voter/Moran dynamics, Eq. ([]]), 
can be coarse-grained by performing an appropriate degree average, to yield 

g{k,k')P{k'\k) 



P w (k -> k') 



E q 9(k,q)P(q\ky 



(3) 



The relevant function defining voter and Moran processes is the probability H(k\ s) that 



18 



-|20|. This 



a spin s at a vertex of degree k flips its value to —s in a microscopic time step 
function can be expressed, within the annealed weighted network approximation, in terms 
of the density Xk of +1 spins in vertices of degree k, taking the form 



IL v (k;+l) = P{k)x k ^P w {k^k'){l 

k> 



Xk', 



(4) 



n v (fc;-l) = P(k)(l-x k )Y,P w (k^k')x k ,, (5) 

k' 
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For 



for the voter model. The origin of these probabilities is easy to understand 
example, Eq. (J4]) gives the probability of flipping a vertex of degree k in the state +1 as the 
product of the probability P(k) of choosing a vertex of degree k, times the probability x k 
that the vertex is in the state +1, times the probability k chooses to interact with a neighbor 
vertex k', which is in state —1 with probability 1 — x k >, averaged over all possible neighbor 
degrees k'. Analogously, the flipping probabilities for the Moran process can be expressed 
as 

n M (A;;+l) = p ( k ')0- - x k >)P w (k' k)x k , (6) 
k' 

Tl M (k; -1) = £ P(k')x k/ P w (k' ->*)(!- (7) 

k' 

Let us now present separately the mean-field analysis for the two models under consid- 
eration. 



A. Voter model 



1. Rate equation, conservation laws and exit probability 



Let us consider the time evolution of the density x k , which is determined in terms of a 
rate equation. Following (5), [lTh^oj], this rate equation is shown to take the form 



Xk(t) = ^P w (k -> k')x k >(t) - x k (t) 
g(k,k')P(k'\k) 

-V^9(k, q )P( q \kf k,it) - Xk ^ (8) 

where in the last expression we have used Eq. ([3]). The complete expression Eq. (jHJ), valid 
for any correlation and weight patterns, is quite difficult to deal with. In order to obtain 
closed analytical expressions, we assume that the underlying network is degree uncorrelated, 



namely, P(k'\k) = k'P(k')/(k) 30|, and moreover, that the weights are simple multiplicative 



functions of the edges' end points, that is, g(k, k') = g s {k)g s {k'). In this way, Eq. (jHJ) becomes 

x k (t) = u v (t) - x k (t), (9) 
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where we have defined 



«*) = E**g3fW do) 



and (/(*)) =E*P(fe)/(fe). 

It is easy to see that the total density of +1 spins, x = J2 k P(k)x k is not a conserved 
quantity, x = — x + uy. The quantity Uy, however, is conserved, uy(t) = 0, as we can see 
by inserting Eq. (j^J) into the time derivative of LOy(t). Finally, the steady state condition of 
Eq. dHJ), x k = 0, implies x k = ojy 

As for the usual voter model [3, Q] the conservation law allows the immediate deter- 
mination of the exit (or "fixation") probability E, i.e. the probability that the final state 
corresponds to all spins in the state +1. In the final state with all +1 spins we have uy = 1, 
while uy = is the other possible final state (all —1 spins). Conservation of uy implies 
then uy = E ■ 1 + [1 — E] ■ 0, hence 

E = UJy. (11) 

Starting from an homogeneous initial condition, with a density x of randomly chosen vertices 



in the state +1, we obtain, since ojy = x, Eh(x) = x as in the standard voter model [ll|. On 
the other hand, with initial conditions consistent of a single +1 spin in a vertex of degree k, 
we have E 1 (k) = kg s {k)/[N(kg s {k))}. 



2. Consensus time 



The backward Fokker-Planck formalism 3l| can be applied to obtain expressions for the 
consensus time TW 

function of the initial density x of +1 spins and the system size 



N. However, followin g R efs. |l7H20l|. it is simpler to apply a one-step calculation and use 



the recursion relation 



18|420] 



T N ({x k }) = Y / T, u v(k; s)[T N (x k - sA k ) + At] + Q({x k })[T N ({x k }) + At], 



:i2) 



where Q = l — J2 s Sfc ^-v(k; s) is the probability than no spin flip takes place, At = 1/N and 
A k = l/[NP(k)] is the change in x k when a spin flips in a vertex of degree k. Rearranging 
the terms in Eq. ffT2l) . and expanding to second order in A*., we obtain 



2N 



E 



X k + UJy — 2x k CUy d 2 T L 



N 



P(k) 



dxl 



dT, 



N 



dx k 



■1. 



(13) 
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Eq. (1131) is simplified by observing that 
in two well distinct temporal regimes 



;he ordering dynamics of the voter model is separated 



li 



321 ] . Over a short time the different densities Xk 



all converge from their initial value to the common value at the steady state Xk = ooy. For 
infinite-size systems, this state survives forever. For finite size N, the system enters instead 
a different regime where the dynamics of densities Xk is enslaved by the fluctuations of uy, 
which performs a slow diffusion until it hits the absorbing values or 1. The consensus 
time is dominated by this second regime. This allows to apply the steady state condition 
which cancels the drift term in Eg . (fT31) . Taking as relevant quantity the conserved weighted 



magnetization ujy, we obtain 
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1 , {k 2 g s {k) 2 ) d 2 T N 



1. 



(14) 



The integration of this equation leads to 

(kg s (k)) 2 



T N {uy) 



-N- 



[uylnuy + (1 — uiy) ln(l — Uy)] 



(k 2 g s (m 

Thus, the ordering time starting from homogeneous initial conditions, Xk = 0Jy 

(kg s (k)) 2 



1/2 is 



T N (x = 1/2) = JV(ln2) 



(k 2 g s (k) 2 ) 



(15) 



B. Moran process 



1. Rate equation, conservation laws and exit probability 



The derivation of the rate equation for the density Xk in the Moran process follows the 
same steps as in the voter model, taking the form 

V k E q g{k\q)P{q\k') 

where in the last step we used the degree detailed balance condition kP(k)P(k'\k) = 
k'P(k')P(k\k') 33J. Assuming again a degree uncorrelated network, and multiplicative 
weights, we are led to 

(16) 



^(t) = P^-[x k (t) - x(t)}. 



(kg s (k)} 
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Again, the total density of +1 spins, x = J2k P(k) x k is n °t conserved, while instead the 
quantity 

1 ^ P(k) 

WM "M1t* (17) 

is conserved, Ojm = 0. Finally, from the steady state condition, x k = 0, we obtain x k = x. 
From the conservation of Um the exit probability is immediately derived as 

E = u M . (18) 

Homogeneous initial conditions lead again to E^ix) = x, while a single +1 spin in a vertex 
of degree k leads to 

El{k) = kg s (k)N([kg s (k)]^Y (19) 
It is interesting to note that in the conserved quantity of the voter model, cuy, each 
density x k is weighted with the product kg s {k) (Eq. (flQj) ). while in the correspondent ujm 



for the Moran process the weight is precise 
noted in the case of unweighted networks 



y the inverse, namely (kg s (k)) 1 (Eq. (fT7|) ). As 



201 ]. intuitively this can be ascribed to the fact 



that in the voter model it is the first selected node that may change its state, while in the 
Moran process it is the second one. Thus in the voter model small-degree nodes change their 
state more often than high-degree nodes, and weighting them with the probability of being 
chosen (kg s (k)) compensates this disparity leading to the conserved quantity toy. Vice versa 
in the Moran process low-degree nodes change their state less often than high degree nodes, 



and the inverse weighting balances this difference 



2. Consensus time 



20|. 



Following the same steps presented for the voter model, and performing the appropriate 
expansion to second order in we obtain the equation 

x - kg s (k) _ dT N 
^(kg s (k)) [Xk X) dx k 

1 . ^ kg s (k) x k + x - 2x k x d 2 T N 
+ 2N^(kg s (k)) P(k) dx\ = ■ 

The steady state condition, x k = x, leads to the cancellation of the drift term. The diffusion 
term is simplified by changing variables with the conserved quantity %, leading to the 
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equation 

1 1 u M (l-u M )d 2 T N 

N ([kg s (k)]-i) (kg s (k)} dcojj ' 1 j 

where we have used the fact that, in the steady state, x = ujm- The solution of equation for 
the consensus time leads now to 

T N (u) M ) = -Nikgsik^dkgsik)}- 1 ) [u) M hx(u) M ) + (1 - u M ) ln(l - u M )\ ■ (21) 

Thus, starting from homogeneous initial conditions, Xk = x = ujm = 1/2, we have 

T N (x = 1/2) = 7V(ln2)(% s (A;))([^ s (A ; )]- 1 ). (22) 

III. NETWORKS WITH POWER-LAW DEGREE DISTRIBUTION AND 
WEIGHT STRENGTHS 

The actual behavior of the exit probability and the consensus time depends, in view of 
the previous calculations, on the topological properties of the network under consideration, 
as well as on the strength of the weights, as given by the function g s (k). In this Section we 
consider explicitly these dependencies for the particular case of networks with a power-law 
degree distribution form, P(k) ~ fc -7 , and a weight strength scaling also as a power of the 
degree, g s (k) = k e . This last selection is reasonable in view of the weight patterns empirically 
observed in real networks Let us focus on the consensus time with homogeneous (x = 0.5) 
initial conditions for the two models considered. 

In the case of the voter model, the ordering time with homogeneous initial conditions 
and weights scaling as a power of k takes the form 

T^(l/2)=iVln(2)|^. (23) 

From this expression, we can obtain different scalings with the network size N, depending 
on the characteristic exponents 7 and 9; we consider only 7 > 2. Using the fact that 
(k a ) ~ const, for a < 7 — 1 and (k a ) ~ /c" +1 ~ 7 for a > 7 — 1, where k c is the upper 
network cutoff, and, in view of the comparison with numerical results for the Uncorrelated 



Configuration Model 34|, considering the scaling k c ~ N 1 ^ 2 for 7 < 3 and k c ~ jV 1 ^ 7 ^ for 
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FIG. 1: (Color online) Phase diagram of the voter model (left) and Moran process (right) on 
weighted scale- free networks. 



7 > 3 



35| . we obtain the following scaling for consensus the time: 



T N (l/2) 



const. 

jy(7-2*-l)/2 
jy2( 7 -fl-2)/(7-l) 

N 



9> 7-2, 7<3 

# > 7 - 2, 7>3 

7-2 > 9 > (7-3)/2, 7 < 3 

7-2 > 9 > ( 7 -3)/2, 7 > 3 

9 < ( T - 3)/2 



(24) 



In Fig. [T] (left) we represent graphically the different scalings of the consensus time T/v in 
the (9, 7) space. 

For the Moran process, the ordering time scales with the network size through the ex- 
pression 

Tjv(1/2) = N\n(2)(k 1+0 }(k^- } 
For 7 > 2, the different possible scalings are as follows: 



(25) 



T N (l/2) 



jy(l+fi)/( 7 -l) 

N, 

jy(2-9-7)/2 



# > 7 - 2, 7<3 

# > 7 - 2, 7>3 
-7 < 6 1 < 7 - 2 
6* < -7, 7 > 3 

# < -7, 7 < 3 



(26) 



Fig. [U (right) depicts the different regimes associated to the behavior of Tv in the (9, 7) 
space. 
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Some comments are now in order. First we notice that all relevant quantities are in fact 
function of the combination k e+1 . This implies that for 9 — — 1 both voter and Moran 
dynamics are predicted to give the same results at the mean-field level, independently of 
the degree distribution. In fact, 8 = — 1 implies that both interacting vertices are extracted 
completely at random (independently of their degree) so that the asymmetry distinguishing 
the voter model from the Moran process vanishes. For other values of 8, on the other 
hand, the effect of weights appears to be completely different for the two dynamics. For the 
voter model, positive values of 8 tend to reduce the consensus time, while 8 < leads to 
increased Tjy, i.e. the dynamics becomes slower. In any case the consensus time is at most 
proportional to the system size N: The dynamics is always relatively fast. Interestingly, the 
HMF analysis predicts the presence of a region (8 > 7 — 2 and 7 > 3) for which the consensus 
time is constant, i.e. the dynamics undergoes an instantaneous ordering process, in contrast 
with what happens in other regions, in which ordered regions of opposite states can coexist 
for very lon g ti mes, reaching consensus only in finite systems and through a large stochastic 



fluctuation 36] . As it will be shown below, this is true only on annealed networks in which 
the quenched disorder imposed by the actual connections in the network is not considered. 
Numerical simulations performed on quenched graphs give different results. 

For the Moran process, on the other hand, T ~ N represents a lower bound for the 
scaling of the consensus time: The dynamics is always rather slow, with an exponent larger 
than 1 for all (7,$). Remarkably, the scaling of T/v turns out to depend symmetrically on 
|0 + 1|: A large positive or a large negative value of 8 + 1 are equally effective in slowing 
down the ordering process. 

IV. COMPARISON WITH NUMERICAL SIMULATIONS 
A. Algorithms 

In order to check the analytical predictions for the voter model and Moran process, we 
have performed numerical simulations of both models on uncorrelated networks generated 



using the Uncorrelated Configuration Model (UCM) 34J. The networks have a degree ex- 



ponent 7, a minimum degree k m = 4 and a maximum degree smaller than or equal to v N 



preventing the generation of correlations for 7 < 3 35]. A weight strength g s (k) = k is 
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imposed by selecting a nearest neighbor j of a vertex % with probability 

p - = r4jfc*' (27) 

where V(i) is the set of nearest neighbors of i. 

Moreover, since HMF equations describe in an exact way dynamics taking place on an- 
nealed networks [ij, we have simulated the voter model and the Moran process also on such 
structures, in order to provide a benchmark of our analytical results. In annealed networks, 
in fact, all links are rewired at each microscopic time step, so that no dynamical correlation 
can build up, and the absence of correlations assumed by mean-field approaches is actually 
implemented. In weighted networks, the probability that a vertex of degree k interacts with 
a vertex of degree k! is given by 

_ k'g s {k')P{k') _ k ll+e P(k') 
P ^ k ^ k )~ {kgs{k)) ~ {k i +e) > ( 28 ) 

where in the last equality we have assumed again that g s {k) = k e . An annealed weighted 
network is thus implemented by choosing as neighbor of any given vertex another vertex of 
degree k, randomly chosen in the network with probability proportional to k 1+e J^. In a 
quenched network, on the other hand, the neighbors of the first node are of course fixed and 
the choice is restricted to them. 



B. Exit probability 

While for homogeneous initial conditions both the voter model and Moran process lead 
to an exit probability equal to the standard voter model, i.e. E(x) = x, invasion initial 
conditions starting from a single +1 spin in a vertex of degree k lead to exit probabilities 
that depend explicitly on the initial degree considered. In particular, we find 

El° tCT {k) ~ k l+ \ Ef OTan {k) ~ k^ 1+e) . (29) 

While for the voter model a single +1 vertex has better chances to invade the system if 
it starts from a high degree vertex, for the Moran process the situation is precisely the 
opposite, a single +1 spin being favored when initially located in the vertices of smallest 
degree. This kind of behavior is actually to be expected from the very definition of the 



models, and has been already reported in unweighted networks 

12 



20] . In fact, a high degree 
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FIG. 2: (Color online) Exit probability E\(k) starting from a single +1 spin in a vertex of degree k, 
for the voter model (full symbols) and the Moran process (empty symbols). Dashed lines represent 
the expected theoretical scaling with k, circles refer to the case = and squares to 9 = 1. Data 
from quenched networks of size iV = 10 3 with 7 = 2.5 (voter model) and 7 = 2.2 (Moran process). 

is beneficial in the voter model since it corresponds to a larger probability of being chosen 
as a partner by a neighbor in search for an opinion to copy, while in the Moran process 
having many neighbors implies a larger probability to be invaded by the opinion at one of 
them. In Fig. [2] we plot the values of the exit probability E\{k) computed from numerical 
simulations. The results fit quite nicely the mean-field predictions in Eq. f|2Ul) : The larger 
the weight intensity, the stronger the impact of high and low degree vertices in the voter 
and Moran processes, respectively. 



C. Consensus time 



In Fig. |3]we check the validity of the scaling behaviors predicted by the HMF treatment 
and sketched in Fig. [TJ In this figure, we plot the scaling of the consensus time T/v as a 
function of N, for different points in the six regions in which the respective phase diagrams 
are divided, compared with the corresponding theoretical mean-field predictions. 

Fig. [3] shows that, overall, the agreement between the scaling predicted by theory and 
numerical data in annealed networks is, as expected, very good. With respect to the results 
for quenched networks, the agreement between HMF theory and simulations is in general 



restricted to small absolute values of 9, as reported for other dynamical processes 



5|. In 
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FIG. 3: (Color online) Scaling with iV for the voter model (left) and Moran process (right) on scale- 
free weighted networks in different regions of the corresponding phase diagrams, Fig. [TJ Squares 
represent data from simulations run on annealed networks, while circles concern quenched graphs. 
Dashed lines represent the theoretical scaling predicted by HMF theory. 
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FIG. 4: (Color online) Consensus time for the voter model (left) and Moran process (right) in slices 
of the phase diagrams at two fixed values of 7 and varying 9, compared with the corresponding 
HMF predictions, Eqs. (|23|) and (|25p . respectively (full lines). Results from simulations performed 
on quenched networks of size N = 3 x 10 3 . 

order to set better limits to the validity of the HMF approximation, in Fig. H] we report the 
numerical values of the consensus time obtained from simulations in quenched networks of 
fixed size N = 3 x 10 3 in slices of the phase diagrams in Fig. [1] performed at two constant 
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FIG. 5: (Color online) Probability that an edge connected to a node of degree k and selected for 
the dynamics is active as a function of k, in the quasi steady state for x = 1/2. The solid line is the 
result of the heterogeneous pair approximation, while symbols are results of numerical simulations. 
Binning has deliberately been avoided to show the large variability of numerical results for larger 
degrees. Data from voter dynamics on quenched networks of size N = 10 5 , with 7 = 3.25 and 

e = 1.5. 

values of 7, one larger and one smaller than 3, and varying 9. These numerical values are 
compared with numerical evaluations of the theoretical predictions in Eqs. ( )23|) and (1251) . 
From Fig. H]we observe that the HMF approximation yields reasonably correct results except 
for large values of 9 (for 7 > 3) or large values of —9 (for 7 < 3). When these errors occur, 
consensus time is underestimated by HMF for the voter model, while it is overestimated 
for the Moran process. At the present stage we are not able to predict a priori when the 
theoretical results fail to describe the behavior of the dynamics taking place on quenched 
networks, but the numerical evidence suggests that the theory works well for values of \9\ of 
the order of those observed in real networks [l|. 

V. HETEROGENEOUS PAIR APPROXIMATION 

There are several possible assumptions in the HMF treatment which could fail when the 
approach breaks down. One is the assumption that the time to reach consensus is dominated 
by the diffusive wandering of the quasi steady state, which is much larger than the time to 
reach such state. More important is however the possibility that the very first hypothesis at 
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the core of mean-field theory, namely that the dynamics of the system can be fully described 



in terms of the densities x k , breaks down 



2n 



32] . This assumption can be violated at several 



different levels. A mild violation occurs when the probability of a node to be in a +1 state 
is correlated with the state of its nearest neighbors. In order to ascertain this possibility, it 
is useful to consider the quantity p k , defined as the probability that an edge connected to a 
node of degree k and selected for the dynamics is active, i.e. it connects nodes in a different 
state. Focusing on the case of the voter model, HMF theory, which explicitly assumes the 
lack of dynamical correlations between the vertices at the ends of any edge, predicts that 
this quantity should be equal to 

pk = ^P w (k ->■ k') [x k (l - xw) + (1 - x k )x k '\ 
k' 

= x k (l - u v ) + (1 - x k )u v (30) 

and hence, for initial homogeneous conditions x k = x = uiy = 1/2, we should have in the 
stationary state = 2x(l — x) = 1/2. 

Fig. shows that, for a case where mean-field is not accurate, this assumption is not 
correct in two respects: Firstly, the value of p k is lower than 1/2 (dashed line), indicating 
that, in fact, correlations build up in the system. Secondly, p k depends on k, implying that 
those correlations depend moreover on the degree of the nodes. 

In order to take into account these degree-dependent dynamical correlations, one needs 
to consider, as relevant dynamical variable, the probability p k<k > that an edge connecting 
a node of degree k with another node of degree k! is active, i.e. the two nodes are in a 
different state. This approach, termed heterogeneous pair approximation (HPA), has been 
introduced and applied to voter models on unweighted networks in Ref. (28] . To determine 
the equation of motion of the quantity p kjk > in the voter model, we observe that this quantity 
is modified if the flipping node has degree k and one of its neighbors has degree k! (or vice 
versa). Let us assume that the flipping (first selected) node has degree k and call k" the 
degree of the copied (second selected) node. It is useful to consider separately the two cases 
where k" ^ k! or k" = k'. 

In the first case the variation Ap k)k / for a single dynamical step (occurring over a time 
At = 1/N) is determined as follows: The probability that a node in state s and degree k 
flips is given by the probability P{k) that the first node selected has degree k times the 
probability <r(s) that it is in state s, times the probability P w (k — >• k") that the second 
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has degree k" multiplied by the probability p k>k "/[2a(s)] that the link connecting the two is 
active. One has then to multiply this quantity by the associated variation of the fraction of 
active links between k and kl . Among the k — 1 other links of the flipping node, the number 
of those connecting to a node of degree k! will be j distributed according to a binomial 
R(j, k — 1) with probability of the single event equal to P(k'\k). In their turn, only n out 
of these j links will be active, with n binomially distributed (B(n,j)) with single event 
probability p ktk r/[2a(s)}. Finally one has to multiply by the variation of p k) y when n out of 
j links go from active to inactive as a consequence of the flipping of the node in k. This is 
given by the variation of the number of active links [(j — n) — n] divided the total number of 
links between nodes of degree k and k', namely NkP{k)P{k'\k). One has then to sum over 
k" ^ k', s, j and n, obtaining 

Aflk,* = P(fc) 5>(a) £ P w (k^k") f 



3 j — 2n 



By performing explicitly the summations (and using J2s V ^ 5 ) = 4/(1 
2x — 1 is the magnetization) the formula becomes 



m 



P w (k k")p k>k > 

k"+k' 



(k-1) 



(31) 

where m = 
(32) 



At Jrfj.. k V 

When k" = k', the value of Ap^^' is similar to Eq. ( 13~T|) with (obviously) P w (k — > k') instead 
of P w (k — > k"), no sum over k", and in the numerator of the last factor j+1 — (n+1) — (n+1) = 
j — 2n — 1, because there are j + 1 links to nodes of degree k, n + 1 of which are active in 
the initial state and inactive in the final. Summing up the two contributions and adding the 
symmetric terms with k and kl swapped, we get 



dp k ,k' k - 1 k' - 1 



dt 



pk- 



\~Pk' — -r, 1 — Pk,k' 



P w (k -+ k') 1 P w (k' k) 1 2p k k-1 2p k , k' - 1 

P(k'\k) k + P(k\k') k< + 1 -m 2 k + 1 - m? k> 

(33) 



When uncorrelated networks are considered, so that 

P w {k k') (k) _ k , e 



P(k'\k) (k 1+e ) 



(34) 



we are led to the final equation 
dp kk i k — 1 k' — 1 



dt 



Pk- 



Pk> 



Pk,k' 



(k) k' e k 



+ 



2p k k-1 , 2p k , k'-l 



(k 1+e ) \k k' ) 1-m 2 k 1 - m 2 k' 



(35) 
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FIG. 6: (Color online) Comparison of the consensus time T/v as a function of N obtained in 
numerical simulations for voter dynamics on quenched networks (empty circles) and the results of 
the numerical evaluation of the HMF, Eq. (|23p (filled circles), and HPA, Eq. (|37|) (empty squares) 
predictions. Data correspond to networks with 7 = 2.75, 9 = 0.5 (top), 7 = 3.25, 9 = 1.5 (center), 
and 7 = 4, 9 = 4 (bottom). 

where, m is the magnetization and, at odds with the case of unweighted networks, the 
definition of p k is now 

Pk = p w( k ~^ k ')pk,k> (36) 

k' 

Solving numerically this equation in the stationary state, it is possible to determine pf, 
which turns out to be in good agreement with numerical simulations, see Fig. Moreover 
it is possible to compute the consensus time Tn, which for the voter model turns out to be 

T N{kl+r (37) 

N 2E k P(k)k 2 ( l + e )p s k (x = 1/2) ' 1 J 

A remarkable agreement between this expression (evaluated numerically) and simulations 

is found even for some cases where HMF theory fails. Notice that no parameter is fitted. 
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Thus, as we can see in Fig. for small values of 9 (9 = 0.5) and small 7, both HMF theory 
and the HPA provide accurate results for the consensus time. Larger values of the weight 
exponent (9 = 1.5) are well represented by HPA, while HMF fails. 

For larger values of 9, however, even the HPA approximation is not sufficient to capture 
the correct behavior of the model. In this regime, a much harsher breakdown of the HMF 
assumptions occurs (5), Q: The state of a node of degree k (or of an edge joining vertices 
of degree k and k') depends not only on the degrees but on the detailed quenched structure 
of the network, much beyond single-node or single-pair features. For example, as 9 — > 00 
j^J each node interacts deterministically with its most connected neighbor. According to 
HMF equations, which describe an annealed scenario, this means that every node will select 
the most connected node(s) in the network. However, in a quenched structure each node 
can choose its partner only among its neighbors, with the result that different portions of 
the network will effectively become independent from the point of view of the dynamics. 
Different regions of the network may therefore order in different states, and in this case the 
final global consensus will never be reached (see also jsj). 



VI. CONCLUSIONS 



We have presented a detailed investigation of the behavior of voter model and Moran 
processes on weighted complex networks. From the analytical point of view we have put 
forward a theoretical framework that allows to deal with generic edge weights. For a spe- 
cific form of the weights we have derived in detail all relevant properties of the dynamical 
processes, such as the exit probability and the scaling of the consensus time as a function 
of the network size. It turns out that the presence of weights has the effect of slowing down 
the Moran process with respect to the unweighted case, while it generally speeds up order- 
ing with voter dynamics. Numerical simulations are in good agreement with the theory for 
small absolute values of 9, while for large \9\ substantial discrepancies show up. An im- 
proved mean-field- like theoretical approach (heterogeneous pair approximation) taking into 
account two-body correlations gives better agreement with numerics. Still in the limit of 
large positive (negative) 9, when the state of a node tends to be deterministically enslaved 
to the state of its neighbor with largest (smallest) degree, the theoretical approaches fail to 
describe in a satisfactory manner the behavior of the system. 
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The positive news is that the mean-field equations describe quite well the dynamics 
observed in real (quenched) networks for weight intensities of the order of the ones observed 
in real- world networks However, the generality of this finding, as well as the intrinsic 
limits of the theory, are in need of a better understanding (see also ^]). A theoretical 
approach able to take into account the detailed quenched structure of weighted networks is 
in order to successfully tackle this problem. 
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